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Abstract 

The spin-1/2 Ising-Heisenberg tetrahedral chain is exactly solved using its local gauge symmetry, which enables one to 
establish a rigorous mapping with the corresponding chain of composite Ising spins tractable within the transfer-matrix 
approach. Exact results derived for spin-spin correlation functions are employed to obtain the frustration temperature, 
at which a product of correlation functions along an elementary triangular plaquette becomes negative and the relevant 
spins experience a spin frustration. In addition, we have exactly calculated a concurrence quantifying thermal entan- 
glement along with a threshold temperature, above which concurrence as a measure of thermal entanglement vanishes. 
It is shown that the frustration and threshold temperature coincide at sufficiently low temperatures, while they exhibit 
a very different behavior in the high-temperature region when tending towards completely different asymptotic limits. 
The threshold temperature additionally shows a notable reentrant behavior when it extends over a narrow temperature 
region above the classical ground state without any quantum correlations. It is demonstrated that the specific heat may 
display temperature dependence with or without an anomalous low-temperature peak for a relatively strong or weak 
Heisenberg interaction, respectively. 
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1. Introduction 

The term frustration was first introduced into mag- 
netism by Toulouse Q]]. According to the original con- 
cept, the frustration is basically a topological property 
and refers to a negative sign of the product of coupling 
constants implying an incapability of spins to satisfy all 
pair spin-spin interactions between them. One of the sim- 
plest frustrated spin systems is the Ising model on a tri- 
angular lattice with the unique antiferromagnetic nearest- 
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neighbor coupling flj]. For many years, frustrated spin 
systems have been widely studied mainly in connection 
with their unconventional low-temperature behavior and 
reentrant phase transitions (see Refs. |3|,|4l] and references 
cited therein). An alternative definition of the spin frustra- 
tion was introduced by dos Santos and Lyra, which have 
proposed a negative sign of the product of pair correlation 
functions as another useful hallmark of the spin frustra- 
tion Q5|] . This latter definition has an advantage in that it 
provides plausible criterion whether or not the spin sys- 
tem is frustrated at finite (non-zero) temperature as well. 

Recently, quantum spin chains have attracted a lot of 
attention as they provide an useful playground for the in- 
vestigation of thermal entanglement and moreover, they 
represent promising candidates for the quantum informa- 
tion processing [6, 7]. It is worth noticing that the entan- 
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glement is one of the most exciting properties of quantum 
systems, which does not have any classical analog and 
is most clearly manifested at zero temperature. However, 
thermal fluctuations gradually destroy the quantum entan- 
glement, which may finally completely disappear above 
a certain temperature called as the threshold temperature 

A theoretical investigation of the Heisenberg model on 
a two-leg ladder is motivated not only from the theoret- 
ical point of view, but also from the experimental view- 
point lfl3l [Till . Unfortunately, the frustrated Heisenberg 
model on a two-leg ladder cannot be rigorously solved 
quite generally and there exist only few exact results for 
some special cases HHO- Owing to this fact, differ- 
ent powerful analytical and numerical approaches have 
been employed in order to examine the magnetic be- 
havior of the frustrated Heisenberg two-leg ladder (see 
Refs. 1 17, 18, 2^, 21, 22] and references therein). 
On the other hand, the ground state of frustrated Ising- 
Heisenberg ladder can be rigorously attained from a pre- 
cise mapping correspondence with the quantum Ising 
chain of composite spins and this model still displays 
typical quantum features including remarkable quantum 
phase transitions as discussed in detail by Verkholyak and 
Strecka 123L 12411 . Another interesting quantum spin model 
represents the Heisenberg tetrahedral chain developed by 
Niggemann et al. B25I1 . whose ground state and thermo- 
dynamic properties can be obtained from the local gauge 
symmetry revealing the hidden Ising symmetry. In the 
present work, we will exactly treat the analogous Ising- 
Heisenberg tetrahedral chain. 

The outline of this Letter is as follows. The tetrahedral 
Ising-Heisenberg chain is introduced in Sec. 2 together 
with its ground-state phase diagram. In Sec. 3 we present 
basic steps of the calculation procedure suitable for evalu- 
ating finite-temperature properties such as the correlation 
functions, the frustration temperature, the threshold tem- 
perature for bipartite entanglement, as well as, the zero- 
field specific heat. Finally, we will draw our conclusions 
and summarize the most important findings in Sec. 4. 

2. The Ising-Heisenberg tetrahedral chain 

Motivated by the comments given in the introduction, 
let us consider the spin-^ Ising-Heisenberg chain of edge- 
sharing tetrahedra as schematically illustrated in Fig.QJa). 




Figure 1 : (a) A schematic representation of the spin- A Ising-Heisenberg 
tetrahedral chain. Thick (blue) lines correspond to the Heisenberg cou- 
pling J xz , while thin (red) lines correspond to the Ising coupling Ji ; (b) 
the tetrahedral chain can also be viewed as a two-leg ladder with equal 
interactions along the legs and diagonals. 



For easy reference, the investigated spin system will be 
hereafter referred to as the Ising-Heisenberg tetrahedral 
chain and it may be alternatively viewed as the particu- 
lar case of the Ising-Heisenberg two-leg ladder B23l 124 ] 
with a Heisenberg intra-rung interaction and a unique 
Ising interaction along the legs and diagonals, respec- 
tively (see Fig. Q2b)). The Hamiltonian of the spin-| 
Ising-Heisenberg tetrahedral chain is given by 



i=\ 

+*(^ + ^ W )(^ + i+*Ui)}. (D 



where & a j (<x? ) denote spatial components a - {x,y, z} 
of the spin-^ Pauli operator at a site labeled by suffix i 
and a (b), respectively (see Fig. [TJ. Note furthermore 
that the interaction terms J x and J z determine the spatially 
anisotropic XXZ Heisenberg coupling, whereas the inter- 
action term J\ represents the Ising coupling. 

To proceed further with calculations, let us introduce 
three spatial components of the total spin angular momen- 
tum S f = o a ai + o"? . for two spins coupled by the Heisen- 
berg interaction. Using the above definition one simply 
obtains the identity (S ") 2 - \ + 2o"" ; cr£\, which allows 
us to express both parts of the Heisenberg interaction in 
terms of two commuting operators Sj = S; • S; and Sr 



a(«>:. ; 



^ W ) = ^-7p?) 2 -7' (2) 
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and 



J m 2 - J i- 



b,i 2 r " 4 
It is quite obvious from Eqs. (O and © that the Heisen- 
berg interaction can be rewritten with the help of total 
spin angular momentum 5, and its zth spatial compo- 
nent Sp which correspond to conserved quantities with 
well defined quantum spin numbers 5,(5, + 1) and S 1 ? = 
—Si, -Si + 1, . . . ,Si, respectively. Besides, the Ising in- 
teraction as the other interaction term entering the Hamil- 
tonian ([TJ can also be trivially expressed by means of the 
spin operator 5|. Bearing all this in mind, the Hamilto- 
nian of the spin-i Ising-Heisenberg tetrahedral chain can 
be decomposed into two parts < H - 'Ho + 'H' ', whereas 
<#o = U2J X + J z ) and 



Jz - Jx 



Up to unimportant constant term 'Ho, the Hamiltonian ([]]) 
of the spin- 1 Ising-Heisenberg tetrahedral chain can be 
rigorously mapped to the Hamiltonian © of the classical 
chain of composite Ising spins, since the quantum spin 
number for the total angular momentum is either S ; = 
or 1. It is noteworthy that the Ising interaction Ji di- 
rectly determines the effective nearest-neighbor interac- 
tion in the equivalent classical chain of composite Ising 
spins, while the spatial anisotropy in the XXZ Heisenberg 
interaction determines the effective single-ion anisotropy 
acting on the composite spins. It should be empha- 
sized, moreover, that the relevant gauge transformation 
directly brings the effective Hamiltonian (|4]i into a diago- 
nal form. Hence, one easily obtains also a diagonal form 
for the Hamiltonian of the spin-| Ising-Heisenberg tetra- 
hedral chain 



<H =^(2J X + J z ) [Si(Si + 1) + S M {S M + 1)] 



•-./•. S':.S': , + — 



,=i 

J 7 -J 



(5) 



It is worthwhile to remark that the transformation be- 
tween the quantum tetrahedral spin chain and the classi- 
cal chain of composite spins is only possible because the 



Hamiltonian (Q} satisfies the local gauge symmetry when- 
(3) ever the Ising interactions along the legs and diagonals are 
identical. If both the Ising interactions would be different, 
the absence of the local gauge symmetry would prevent 
a rigorous mapping onto the classical chain of compos- 
ite spins even though the investigated spin system could 
still be mapped onto the quantum Ising chain of composite 
spins [23, 24]. However, the latter quantum chain of com- 
posite spins can exactly be treated only at zero tempera- 
ture unlike the former classical chain of composite spins 
(see Sec. 3). It should be also noted here that the Heisen- 
berg tetrahedral chain under the local gauge symmetry has 
been studied by Niggeman et al. B25I1 when combining 
the fragmentation approach with numerical calculations, 
while the relevant high-temperature expansion has been 
developed in Ref. ||2£ 



^ 2.1. Ground- state phase diagram 



The diagonalized form of the Hamiltonian (|5]l can be 
straightforwardly used in order to obtain all possible 
ground states of the spin-^ Ising-Heisenberg tetrahedral 
chain. The Hamiltonian © implies an existence of four 
different ground states, namely, the ferromagnetic state 
(FM), the superantiferromagnetic state (SA), the triplet- 
dimer state (TD), and the singlet-dimer state (SD), which 
are unambiguously given by the following eigenvectors 



iFM>=f]i:> ; , 

,■=1 

N/2 

|SA> = p]i:> 2i _ 1 ®|:> 2i> 
|TD>=n-^(|i> i + |;> i ) ! 



N 



|SD>=fTi(|!> i -|;> i ). 
*=i V2 



(6) 
(7) 
(8) 
(9) 



±i and the 



In above, the upper ± sign applies for cr^ . 
lower + sign stands for erf . = ±1. The corresponding 
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JJ\J,\ 

Figure 2: The ground-state phase diagram in the J X /\J\\ - J z l\Jl \ plane. 
An upper sector of the phase diagram contains FM (SA) ground state for 
the particular case with J\ > (Ji < 0). 



3. Thermodynamics 

In this part, let us calculate finite-temperature proper- 
ties of the spin-| Ising-Heisenberg tetrahedral chain. To 
gain the partition function, it is quite convenient to adopt 
a trace invariance and to calculate the partition function 
using the diagonalized form of the Hamiltonian (O 



Zn = Tre 



-f¥H 



Tre 



(ID 



Here, j3 = l/ik^T), is being the Boltzmann's constant 
and T is the absolute temperature. Let us first rewrite the 
diagonalized form of the Hamiltonian (0 into the follow- 
ing compact form 



energy per unit cell for each individual ground state reads 



E FM 






E SA 




J z 

J + J 


E TD 


Jz 

~T 


J x 
2 ' 


E SD 


~~4 


+ — . 

2 



(10) 



By making use of the eigenenergies ( TTOb we have con- 
structed the zero-temperature phase diagram, which is 
displayed in Fig. [2]in the J X I\J\\— JJ\J\\ plane. The disor- 
dered SD state with singlet dimers placed on all Heisen- 
berg bonds becomes the ground state whenever J x < 
and J z < —J x — 2\J\\. On the other hand, the similar dis- 
ordered TD state is being the respective ground state if 
J x > and J z < J x - 2\Ji\. The most fundamental dif- 
ference between TD and SD states lies in the symmetric 
versus antisymmetric quantum superposition of two an- 
tiferromagnetic states as evidenced by the eigenvectors 
([8]l and (O, respectively. Note furthermore that TD and 
SD ground states are apparently of a purely quantum na- 
ture and they arise from a mutual competition between the 
Heisenberg and Ising interactions. Finally, the less inter- 
esting classical FM (SA) state without any quantum cor- 
relations constitutes the ground state for Ji > (J\ < 0) 
on assumption that J z > \ J X \ - 2\J\\. 



(12) 



i=l 

with the Hamiltonian defined as 
J, - J 



n 



-■^[Si(S( + l) + S w (S w + l)]. (13) 



Then, the partition function can be partially factorized if 
substituting Eq. ( fl2l i into ( ITTb 



T ,n< 



-m 



i=l 

(14) 

and the latter expression behind the product symbol can 
be subsequently identified as the usual transfer matrix 



2 2 
3 

x y 
xy 

2 2-1 

V xry z 



3 2 2 -1 

x y xy x y z 

4 2 3 

x x x y 

x 2 I xy 

3 2 2 

x y xy x y z 



(15) 



with x = eP J *l 4 , y = eP J --l A , and z = eP J ' . Hence, it follows 
that the partition function can be evaluated by the standard 
transfer-matrix approach M26I1 . Moreover, the eigenvalue 
problem for the transfer matrix (15[ can easily be solved 
from the secular determinant, which leads to the following 
quartic equation 

A (A 2 - a x A + fl ) [A - x 2 y 2 (z - z~ 1 )] =0 (16) 



4 



with 



a Q =x Z y Z z~\z- 1)V + 1), 
a\ -x A + 1 + x 2 y 2 (z + z~ l ). 



(17) 
(18) 



After a straightforward calculation one finds all four 
transfer-matrix eigenvalues 



A\ =\ {ai - ^ - 4a j , 

/b =0, 



whereas the largest eigenvalue is /to by inspection. In the 
thermodynamic limit, the free energy per unit cell is sim- 
ply given by the largest eigenvalue of the transfer matrix 



f = -p- 1 lim T>Zn = 

N—>cc jy 



2J X + J z 1 



-^in^o. (19) 



It is worth mentioning that the first constant term of the 
free energy was obtained during the local gauge transfor- 
mation and this term is irrelevant for most thermodynamic 
quantities. After a straightforward but a little bit cumber- 
some manipulation, one may derive the following final 
formula for the free energy per unit cell 



-El lfiJx\ 
e 4 cosh(j67i) + e 4 cosh y—^-J + Q 



f = -/T' In 
where the coefficient Q is defined as follows 



(20) 



Q 



e ^ cosh(/?7i) 



+ 4 cosh 



m 

(21) 



3.1. Correlation functions 

Now, let us calculate a few pair correlation functions, 
which might be helpful in characterizing the magnetic be- 
havior at finite temperatures. The short-range correlation 
functions can in turn be calculated by differentiating the 
free energy (l20l > with respect to the relevant coupling con- 
stant. For instance, both spatial components of the corre- 
lation function between two spins coupled by the Heisen- 
berg interaction can be obtained by differentiating the free 




Figure 3: Thermal variations of the pair correlation functions for the 
isotropic antiferromagnetic Heisenberg interaction J x = J z < 0, the 
antiferromagnetic Ising interaction J\ < 0, and two different values of 
the interaction ratio: (a) JJ\J\\ = -0.9; (b) JJ\J\\ = -1.1. 



energy < f20b with respect to the relevant spatial component 
(J z or J x ) of the Heisenberg interaction 



(cr 2 cri ) = 
x a,i b.r 



1 

4G 



El _&Z lPJx\ 

e 4 cosh(yS7i) - e 4 cosh^-^-j 



(22) 



and 



/ a-x - x \ i &y a v \ 



sinh(^) [2 - cosh^) + e-^cosh(^) + e J Tg] 
W,r -T cosh(j67i) + e~T i cosh(^) + Q\ 

(23) 



4G 



Due to the local gauge symmetry of the tetrahedral chain, 
the longitudinal pair correlation functions between the 
nearest-neighbor spins interacting via the Ising interac- 
tion must be equal to each other. Thus, the differentiation 
of the free energy (1201) with respect to the Ising interaction 
J\ gives us another important pair correlation functions 

<^,,^,,+i> = <^aL + i> = ^ b ,iK, i+ i) = <^,i^, i+ i> 

e-r sinhOe/i) [e^ coMPJ\) - e-^cosh(^) + Q] 



4g[e^ cosh(y37i) + e-T i cosh(^) + g] 



(24) 



Let us analyze in detail typical temperature depen- 
dences of all calculated pair correlation functions for the 
most interesting particular case with the antiferromag- 
netic Ising and Heisenberg interactions. For illustration, 
the correlation functions are plotted in Fig.[3]against tem- 
perature for the isotropic Heisenberg interaction and two 
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different values of the interaction ratio JJ\J\\ = -0.9 
and -1.1, respectively. Two selected values of the in- 
teraction ratio J Z I\J\\ has been chosen so as to achieve 
two different ground states. As a matter of fact, the 
zero-temperature limits of the displayed correlation func- 
tions are consistent with a presence of the classical SA 
ground state in the former case with J Z J\J\ \ = -0.9, while 
the existence of the more striking quantum SD ground 
state is evident in the latter case with J Z I\J\\ = -1.1. It 
is quite obvious from Fig. [3ja) that thermal fluctuations 
may induce quantum correlations above the classical SA 
ground state if the interaction ratio JJ\J\\ is selected suf- 
ficiently close to the ground-state boundary with the SD 
state. It actually turns out that the transverse pair corre- 
lation function (&* trf ) becomes non-zero upon increas- 
ing temperature and the longitudinal pair correlation func- 
tion (cF .cr| .) changes its signal at some finite tempera- 
ture Tf called as the frustration temperature. Following 
the definition by dos Santos and Lyra Q5D, the tetrahedral 
chain is frustrated above the frustration temperature as 
the product (cr cf, XcF cr . , VcH ■ ,cn ) between three 

* N a,i b,i ' x a,i a,(+l b,v 

correlation functions along the triangular face of an el- 
ementary tetrahedron becomes negative. It is worthy of 
notice, moreover, that the correlation functions {&~ a fi*. ; ) 
and {& z a fr z a remain negative for arbitrary temperature 
whenever the SD state constitutes the ground state. This 
result is taken to mean that the frustration persists for all 
temperatures if the SD state constitutes the ground state. 

3.2. Thermal entanglement 

Another interesting aspect of the spin-i Ising - 
Heisenberg tetrahedral chain surely represents a study of 
thermal entanglement with regard to its unusual quantum 
ground states. For this purpose, one may take advantage 
of the quantity referred to as the concurrence II271 12811 . 
which may serve as a measure of the bipartite entangle- 
ment. The concurrence is defined in terms of the reduced 
density matrix p 

C(p) = max{0, 2A max - TiR], (25) 

where X max is the largest eigenvalue of the matrix 

r = ^Ip&y ®&yp*&y ®&y. (26) 

Here, p* is the complex conjugate of the reduced density 
matrix p and & y is being the relevant Pauli matrix. It is 



quite well established that matrix elements of the reduced 
density operator can be expressed through the correla- 
tion functions 1 28 , 29| and hence, the concurrence can be 
directly calculated from the correlation functions as evi- 
denced by Amico and co-workers j2^l . In our case, the 
concurrence quantifying the bipartite entanglement be- 
tween two spins coupled by the Heisenberg interaction 
then follows from the relation 



C = max -|0,4 \{d*^)\ - 



- + 2<<x z .<xr > 

9 x a,i b,v 



(27) 



which just involves two pair correlation functions given 
by Eqs. (1221 and (|23| >. It is noteworthy that an equiva- 
lent result could be obtained by employing the approach 
developed in Ref . 13011 . 

The concurrence is plotted against temperature in 
Fig. Ufa) for a specific choice of the antiferromagnetic 
Ising and Heisenberg interactions, which ensures occur- 
rence of the classical SA spin alignment in the ground 
state. It surprisingly turns out that the Ising-Heisenberg 
tetrahedral chain exhibits a weak thermal entanglement in 
spite of the classical SA ground state whenever a relative 
strength between the Heisenberg and Ising couplings is 
sufficiently close to the phase boundary Jz/\Ji\ with the 
other possible SD ground state {j\j\J\\ = -1.0 if assum- 
ing the isotropic Heisenberg interaction J x = J z ). Accord- 
ingly, the concurrence shows a striking reentrant behavior 
when it becomes non-zero at a lower threshold temper- 
ature and disappears at an upper threshold temperature 
k&TJ\J\\ ~ 0.4, whereas the lower threshold tempera- 
ture gradually tends towards zero as the interaction ratio 
approaches the boundary value J h /\J[\. It can be readily 
understood that the observed weak thermal entanglement 
comes from a temperature-induced strengthening of the 
transverse correlation function (fr x cr? .), whose effect is 

x a.i b.i' 7 

also supported by a relevant crossover of the longitudi- 
nal correlation function (<r J .en .) from the ferromagnetic 
regime to antiferromagnetic one after passing through the 
frustration temperature. 

In addition, Fig.Hfb) displays more typical temperature 
variations of the concurrence on assumption that the inter- 
action ratio between the Heisenberg and Ising interactions 
supports the SD ground state. The spins coupled by the 
Heisenberg interaction are consequently maximally en- 
tangled at zero temperature and their quantum entangle- 
ment is gradually suppressed by thermal fluctuations. It 
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(a) K T/ \ J ,\ M k n T, \ J , 



Figure 4: Temperature dependences of the concurrence for the isotropic 
antiferromagnetic Heisenberg interaction J x = J z < 0, the antiferromag- 
netic Ising interaction J\ < 0, and several values of the interaction ratio: 
(a)7 z /|/i|>-1.0; (b)/ z /Ui|<-1.0. 

appears worthwhile to remark that the threshold tempera- 
ture above which the relevant spins become disentangled 
becomes higher as the magnitude of the ratio J Z /\J\\ be- 
tween the Heisenberg and Ising interactions increases. 

3.3. Threshold versus frustration temperature 

Now, it might be quite interesting to compare the 
threshold temperature, delimiting the regimes of entan- 
gled spins (finite concurrence) and unentangled spins 
(vanishing concurrence), with the frustration tempera- 
ture that delimits the regimes with (non frustrated) and 
without (frustrated) a well defined spin configuration 
among the three spins forming a triangular face of an 
elementary tetrahedron. It is quite obvious from our 
previous analysis of the correlation functions that the 
absence of spin correlations can be directly connected 
with a disappearance of the longitudinal pair correla- 
tion function {d* p\ •) between two spins coupled by the 
Heisenberg interaction. While the threshold temperature 
can be calculated from the condition C — by mak- 
ing use of Eq. d27l for the concurrence, the frustration 
temperature can be either calculated from the condition 

(K&jXKKm^m^? = or the sim P ler ec l uiv - 
alent condition (o^.oi.) = after taking into account 

Eq. (l22b for the relevant correlation function. 

The frustration and threshold temperatures are plotted 
in Fig. [5] against the interaction ratio J Z /\J\\ for three dif- 
ferent values of the exchange anisotropy J X /\J Z \ of the an- 
tiferromagnetic Heisenberg interaction. It is worthwhile 
to remark that the Ising-Heisenberg tetrahedral chain is 
thermally entangled just in the parameter region below the 
line of threshold temperatures (solid lines), which shows a 
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Figure 5: The threshold temperature (solid lines) and frustration temper- 
ature (broken lines) as a function of the interaction ratio J Z I\J\ \ for the 
antiferromagnetic Ising (J[ < 0) interaction and three different values of 
the exchange anisotropy J X /\J Z \ = -0.5, -1.0 and -2.0. 

relatively narrow reentrant behavior in a close vicinity of 
the ground-state boundary between the SD and S A ground 
states in agreement with our previous analysis of the con- 
currence. On the other hand, the Ising-Heisenberg tetra- 
hedral chain is affected by a spin frustration only above 
the line of frustration temperatures (broken lines), which 
also starts from the ground-state boundary between the 
SD and SA ground states but asymptotically tends to- 
wards infinity as J Z /\J\\ — > 0. Altogether, the frustration 
and threshold temperature remarkably coincide at suffi- 
ciently low temperatures, while they exhibit a very differ- 
ent behavior in the high-temperature region when tending 
towards completely different asymptotic limits. It is quite 
interesting to notice, moreover, that the thermally entan- 
gled region above the quantum SD ground state as well 
as the more peculiar thermally entangled region above 
the classical SA ground state fall into the frustrated re- 
gion. Hence, it could be concluded that the spin frustra- 
tion seems to be an essential ingredient, which is needed 
in order to produce the thermal entanglement in the Ising- 
Heisenberg tetrahedral chain. 

3.4. Specific heat 

Last but not least, let us turn our attention to a dis- 
cussion of typical temperature dependences of the spe- 
cific heat numerically calculated from the free energy < f20b 
with the help of thermodynamic relation C = —NT^. 
Fig. |6ja) displays typical thermal variations of the spe- 
cific heat for the isotropic antiferromagnetic Heisenberg 
coupling J x — J z < by choosing the interaction ratio 



7 



1.5 
1.2 

J 0.9 

O 

0.6 
0.3 



0.0 
0.0 

(a) 





JJ\J, 


= -0.9 




— //IJ, 


= -1.0. 




A //[J, 











0.2 



0.4 0.6 



0.8 1.0 




Figure 6: (a) Thermal vai'iations of the specific heat for the isotropic an- 
tiferromagnetic Heisenberg interaction J x = J, < by selecting the in- 
teraction ratio J Z I\J\ I at and close to the ground-state boundary between 
the SA and SD ground states; (b) Thermal variations of the specific heat 
for the fixed value of the interaction ratio J-I\J\ | = -3 and three different 
values of the exchange anisotropy J X /\J Z \ m the Heisenberg interaction. 



J Z /\J\\ close to and precisely at the ground-state bound- 
ary between the SA and SD ground states. Under this 
condition, the specific heat generally exhibits the tem- 
perature dependence with a single round Schottky-type 
maximum irrespective of a relative strength JJ\J\\ be- 
tween both the interaction parameters, whereas there does 
not appear any additional low-temperature maximum that 
would reflect thermal excitations between the ground state 
and the respective low-lying excited state (either from SA 
to SD state or vice versa). Contrary to this, the more 
diverse temperature variations of the heat capacity can 
be observed by considering the more general case of the 
anisotropic Heisenberg coupling (see Fig.|6fb)). In addi- 
tion to the standard thermal dependences of the specific 
heat with a single more or less symmetric round maxi- 
mum, one may also detect a more pronounced double- 
peak temperature dependence as the exchange anisotropy 
strengthens (i.e. J X /\J Z \ — > 0). It should be pointed out 
that the Ising-like exchange anisotropy drives the investi- 
gated spin system close to the ground-state boundary be- 
tween two quantum SD and TD ground states and thus, 
the additional low-temperature peak clearly reflects the 
low-lying thermal excitation from the SD ground state to- 
wards the excited TD state. 



4. Conclusion 

In this work, the spin-| Ising-Heisenberg tetrahedral 
chain has been exactly solved by employing its local 
gauge symmetry and the transfer-matrix technique. It is 



worthwhile to recall that the investigated spin system can 
be alternatively viewed as the spin-^ Ising-Heisenberg 
two-leg ladder with equal Ising interactions along the legs 
and diagonals, respectively. The elaborated rigorous pro- 
cedure enabled us to obtain exact results for basic ther- 
modynamic quantities (e.g. free energy, specific heat) 
and several pair correlation functions, which were subse- 
quently used for a calculation of the concurrence quanti- 
fying the thermal entanglement between two spins cou- 
pled by the Heisenberg interaction, the frustration and 
threshold temperature. It has been demonstrated that 
the frustration and threshold temperature coincide at low 
enough temperatures, while they tend towards completely 
different asymptotic limits in the high-temperature range. 
It actually turns out that the frustration temperature di- 
verges as the interaction ratio J z /\Ji\ — > 0, while the 
threshold temperature goes to infinity if J z /\Ji\ — > -co. 
An interesting observation is that the investigated Ising- 
Heisenberg tetrahedral chain is thermally entangled just 
in the frustrated region and hence, the frustration seems to 
be an essential ingredient in order to produce the quantum 
entanglement. It is worth noticing that the threshold tem- 
perature for the pairwise entanglement exhibits a striking 
reentrant behavior. Therefore, it can emerge at moder- 
ate temperatures above the classical ground state without 
any quantum correlations. This feature unveils that the 
usual view point of thermal fluctuations as a mechanism 
through it quantum correlations are destroyed is not gen- 
erally correct. In systems with an ordered classical ground 
state but close to a frustration point, thermal fluctuations 
can firstly destroy the magnetic local order, thus leaving 
space to the emergence of quantum entanglement at fi- 
nite temperatures. Only at sufficiently high temperatures, 
the absence of quantum entanglement is recovered. Be- 
sides, we have also investigated in particular typical tem- 
perature variations of the specific heat, which may exhibit 
a thermal dependence with (without) an anomalous low- 
temperature peak for a relatively strong (weak) Heisen- 
berg interaction. Finally, it should be also mentioned that 
our rigorous procedure can be adapted in a rather straight- 
forward way to account for the non-zero external mag- 
netic field as well. Our future work will continue in this 
direction. 
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